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ABSTRACT 



Aims. We explore the morphology of hard (18-30 keV) X-ray emission from the Coma cluster of galaxies. 
Methods. We analyse a deep (1.1 Ms) observation of the Coma cluster with the ISGRI imager on board the INTEGRAL satellite. 
Results. We show that the source extension in the North-East to South-West (SW) direction (~ 17') significantly exceeds the size of 
the point spread function of ISGRI, and that the centroid of the image of the source in the 18-30 keV band is displaced in the SW 
direction compared to the centroid in the 1-10 keV band. To test the nature of the SW extension we fit the data assuming different 
models of source morphology. The best fit is achieved with a diffuse source of elliptical shape, although an acceptable fit can be 
achieved assuming an additional point source SW of the cluster core. In the case of an elliptical source, the direction of extension of 
the source coincides with the direction toward the subcluster falling onto the Coma cluster If the SW excess is due to the presence of 
a point source with a hard spectrum, we show that there is no obvious X-ray counterpart for this additional source, and that the closest 
X-ray source is the quasar EXO 1256+281, which is located 6.1' from the centroid of the excess. 

Conclusions. The observed morphology of the hard X-ray emission clarifies the nature of the hard X-ray "excess" emission from the 
Coma cluster, which is due to the presence of an extended hard X-ray source SW of the cluster core. 

Key words. Galaxies: clusters: Coma Cluster - X-rays: galaxies: clusters - Gamma rays: observations 



1. Introduction 

Clusters of galaxies are the biggest bound structures of the uni- 
verse, and, according to the hierarchical scenario of structure for- 
mation, the latest to form. They are filled by a hot (10^ - lO' K) 
plasma, called Intra Cluster Medium (ICM), and thus radiate in 
soft X-ray bands through thermal Bremsstrahlung. 

Since clusters of galaxies are the latest and biggest struc- 
tures to form, we expect some of them to be still forming, and 
experiencing major merging events with smaller clusters. This 
is the case of the Coma cluster, that is currently merging with 
the NGC 4839 group. In such events, the merging of the ICM 
of the two clusters creates shock fronts in which theory pre- 
dicts that a large population of particles would be accelerated 
to high energies (ISarazinlll999l) . This phenomenon should then 
produce a reheating of the gas and create a higher tempera- 
ture plasma that would radiate more strongly in hard X-rays. 
Alternatively, interaction of the population of mildly relativis- 
tic electrons that produce t he halos of galaxy clusters via syn- 
chrotron radiation (Feretti & Giovanninil 2007h with the Cosmic 
Microwave Background would then produce hard X-ray emis- 
sion through inverse Compton processes, and thus add a power- 
law tail to the spectrum in the hard X-ray domain. Another 
possible model involves a population of multi-TeV electrons 
that would radiate in hard X-rays through synchrotron emis- 
sion (iTimokhin et al. 2004). Detection of this hard X-ray excess 
would help in learning more about the cosmic ray population 
detected by radio observations. Furthermore, characterization of 
the morphology of the hard X-ray emission would bring a pos- 



sible identification of acceleration sites, and since clusters of 
galaxies are one of the few possible candidates for acceleration 
of cosmic rays at high energies, it would bring important infor- 
mation on the origin of cosmic rays. 

Recent reports of detection of a hard X-ray excess 
by Beppo-SAX jFusco-Femiano etaH |2004|) and RXTE 



jRephaeli & Grubeil 120021) 
confirm the existence of a 



in the Coma cluster appear to 
high energy tail of the spectrum 
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of merging clusters, and thus prove the existence of particle 
acceleration sites in these clusters. However, these detection s 
are quite weak and controversial dRossetti & Molendil |2004|) . 
and since the hard X-ray instruments on both Beppo-SAX and 
RXTE are non-imaging, contamination by very hard point 
sources inside the cluster cannot be excluded (e.g. by the central 
galaxy NGC 4874, NGC 4889 or the QSO EXO 1256-H281). 
In addition, no information on the morph o logy o f the hard 
X-ray emission was obtained. iRenaud et alj (l2006al) presented 
an analysis of a first 500 ks set of INTEGRAL data and were 
not able to confirm the presence of a hard X-ray excess or to 
sensibly constrain the hard X-ray morphology of the source. 

In this paper, we use the imaging capabilities of the 
IBIS/ISGRI instrument dLebrun et al.ll2003h to extract informa- 
tion on the hard X-ray emission of the Coma cluster. In Sect. 
|3] we present the results of our imaging analysis of the ISGRI 
data, and compare them with existing XMM-Newton data in the 
soft X-ray domain. In Sect. 14. 11 we describe a method to analyse 
extended sources with a coded mask instrument to extract quan- 
titative flux measurements, and apply it to the case of the Coma 
cluster. In Sect. |5] we present a combined XMM/INTEGRAL 
spectrum of the cluster. Finally, discussion of our results is pre- 
sented in Sect.|6] 
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2. Data 

Our analysis covered 401 Science Windows (ScWs) of public 
data, for a total of 1 . 1 Ms of observation. We analysed the data 
with the latest release of the Offline Scientific Analysis (OSA), 
version 6.0, and eliminated ScWs with a high background level. 
We used the remaining data to create a mosaic image in the stan- 
dard way. Table[T]gives the log of the observation. 



INTEGRAL 
revolution number 


Observation 
dates 


No. of 
pointings 


Observing 
time [ks] 


0036 


Jan 29-31, 2003 


63 


140.1 


0071-72 


May 14-18, 2003 


135 


304.5 


0274-75 


Jan 10-15, 2005 


57 


202.4 


0317-18 


May 19-25, 2005 


99 


333.4 


0324-25 


Jun 9-11,2005 


47 


164.5 




1,144.9 



Table 1. INTEGRAL observation log on the Coma cluster 



The XMM-Newton image is produced using the data of the 
PN car nera taken in Jun e 2000 during the Coma mosaic obser- 
vation jBriel et al.ll2001 ). We used the SAS software version 6.5 
and th e background substract ion method from the Birmingham 
group jRead & Ponmanll2003h to analyse the data. 




Fig. 1. Standard OSA 6.0 significance image of the Coma cluster 
with Ms of data in the 18-30 keV energy band. Significance 
contours from 3 to lOcr in steps of Icr are overlayed in red. The 
position of the 3 brightest X-ray point sources is shown. For 
comparison, the inset in the bottom right corner shows a mosaic 
image of a known point source in the same field, NGC 4388. 



3. Imaging analysis in the hard X-ray domain 

3.1. Mosaic image of tlie Coma cluster 

Figure [1] shows a mosaic image in the 18-30 keV band extracted 
with the standard OSA 6.0 tools, using the data described in 
Table [T] with 3-10 significance contours. For comparison, we 
show in the inset in the right bottom corner of the image a mo- 
saic image of a point source (NGC 4388) produced using a com- 
parable amount of data and normalized so that the amplitude of 
the brightest pixel is the same as in the case of Coma cluster. 

A first look at the image indicates that the Coma cluster 
source is extended. It is not surprizing, since the angular size 
of the cor e of the of Coma cluster in 1-10 keV energy band is 
D ^ 20' jSchuecker et al.| [2004). which is larger than the size 
of the point spread function (PSF) of ISGRI (2cr = 12', which 
is the angle subtended by a 11.2 mm mask hole at a distance of 
3200 mm). 

A comparison of the 18-30 keV image with the image in 1- 
10 keV band obtained with XMM-Newton is shown in Fig. |2] 
The left panel of this figure shows the INTEGRAL image with 
the overlayed contours obtained by smoothing the XMM-Newton 
image with the ISGRI PSF (modelled as a gaussian of full width 
12'). The right panel of the figure shows the residuals of the 
INTEGRAL mosaic image after substraction of the smoothed 
XMM-Newton image, renormalized in a way that the difference 
between INTEGRAL and XMM-Newton flux cancels at the max- 
imum of the XMM-Newton emission. One can clearly see that 
significant residuals are left in the South- West (SW) part of the 
INTEGRAL source after the substraction. This indicates that the 
hard X-ray source detected by INTEGRAL is more extended in 
the SW direction than the XMM-Newton source. 



3.2. Source morpliology from mosaic image 

To clarify the nature of the SW extension of the INTEGRAL 
source we attempted to fit the ISGRI image assuming different 



surface brightness distributions /(r) of the hard X-ray source. 
Specifically, we consider the following possibilities: 

- Model 1: A single point source given by a Gaussian with the 
half-width equal to the size of the PSF of ISGRI, 



/(r) = A exp 



2o-2 



(1) 



where ro is the central vector position of the source which is 
left free while fitting; 

Model 2: A superposition of two point sources with overlap- 
ping PSFs, 



/(r) = A exp 



(r-ri)2 



2cr2 



-I- Bexp 



ir-r2f 



(2) 



where ri , ri are the positions of the two sources which are 
left free; 

Model 3: An ellipse-shaped extended source with the sur- 
face brightness profile 



I{r) - A exp 



((r - ra) • nf ((r - x nf 



20-? 



2^2 



(3) 



where n is the unit vector in the direction of the major axis 
of the ellipse, is the position of the centroid of the ellipse 
and (Ti , cTi are the sizes of the major and minor axes of the 
ellipse which are all left free while fitting; 
Model 4: A superposition of an extended source with mor- 
phology of the core of Coma cluster in the 1-10 keV energy 
band (the surface brightness profile described by a Gaussian 
convolved with the ISGRI PSF, i.e. a Gaussian with a half- 
width of 10') and of an additional point source 



I(r) - A exp 



2cr\ 



+ B exp 



2cr2 



(4) 
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Fig. 2. Left: INTEGRAL mosaic image with contours from XMM-Newton overlayed. The XMM-Newton image is smoothed with a 
gaussian of the width 12' to match the angular resolution of INTEGRAL. Right: the residuals after the subtraction of the XMM- 
Newton profile from the INTEGRAL image (see text). The South- West excess in the residual image is apparent. 



where rc is the position of the centroid of the soft X-ray emis- 
sion and 0-3 - 10'. r4 is the position of the additional point 
source which is left free while fitting, and cr is fixed to the 
half-width of the PSF in the same way as in model 1 . 

We fitted the 21 x21 pixels (105'x 105') part of the image around 
the catalog position of the Coma cluster minimizing the;t'^ of the 
fit, defined as 

2 ^ y (IMG(/,;)-/(/,;))^ 
^ WARdJ) • 

where IMG(!,j) and VAR(/, 7) are the values of intensity and 
variance in a given image pixel (/, j). The best fit results for the 
four models are shown in Fig. [3] 

Fitting the intensity image with a single point source (Model 
1), we find that the source position (the vector ro in ([T])) is shifted 
compared to the centroid of the XMM-Newton image in the SW 
direction by ARA= 3.9', ADEC= 0.5'. The rather high reduced 

of the fit with a point source model, x^^d ~ ^-^^ ^^'^^ '^■^^ 
degrees of freedom) indicates that the single point source model 
does not provide a good description of the source morphology, 
confirming the analysis shown in Fig. 2. 

Fitting the source morphology with a two point sources 
(Model 2), one finds a better reduced x^^.^ - 1-33. The two 
point sources model provides the possibility of finding the di- 
rection of extension of the source. Namely, the best fit is pro- 
vided by the model in which the two point sources are sit- 
uated at RAi = 195.04 + 0.01, DECi = 28.00 + 0.01 and 
RA2 = 194.77+0.01, DEC2 = 27.90±0.01. The angular distance 
between the two point sources is di2 - \ri - r2\ = 15.3' ± 0.6'. 
The ratio of intensities of the two sources is A/B = 0.84. The 
distance between the two point sources is larger than the size of 
the PSF of ISGRl, which confirms again that the source cannot 
be described by Model 1 . 

The best fit to the morphology of the source is found when 
fitting the image with the model of an elliptically-shaped source 




Fig. 3. Comparison between the best fits to INTEGRAL mosaic 
image (see Fig. [U with a standard point source PSF (left), the 
PSF of 2 point sources (second from the left), the PSF of an 
extended source with the shape of an ellipse (third from the left) 
and the PSF of a source with the morphology of the Coma cluster 
in the 1-10 keV plus an additional point source (right). 

(Model 3). The fit results in a reduced = 1.23. The pa- 
rameters of the best fit model are cri - 16.8' + 0.5', 0-2 = 
11.7' + 0.4' and the coordinates of the centroid of the ellipse, 
r3,RA=194.89 + 0.01, DEC= 27 .94 + 0.0 1 . The direction of the 
major axis of the ellipse, n, is inclined at the angle 6 = 61 +4°. 
One can see that the fitted position of the centroid of the ellipse 
is shifted in the same direciton as the position of the single point 
source fit, ro. The size of the major axis of the ellipse, cr; is 
roughly equal to dn- 

The fit by an extended source of the shape of the Coma clus- 
ter in the 1-10 keV band plus an additional point source (Model 
4) gives a good fit, x^^j - 1-36. The additional point source is 
found at RA=194.7l'+ 0.01 and DEC=27.87 + 0.01, which is 
located 6.1' away from the quasar EXO 1256H-281. Fig.|4]shows 
the position of the fitted point source on the residual image with 
Icr, 2o- and 3o- error contours. 

4. A method to analyse extended sources with a 
coded mask instrument 

4. 1 . Description of the method 

The imaging capabilities of ISGRI make it possible to obtain 
information on the morphology of the hard X-ray emission of 
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Fig. 4. Residual image after substraction of the scaled XMM- 
Newton image from the INTEGRAL mosaic image. The position 
of the fitted point source from Model 4 is displayed with Icr, 
2cr and 3cr error contours. The positions of the brightest AGNs 
emebedded in the cluster is also displayed. 



the Coma cluster and other slightly extended sources for the first 
time. However, the standard Of fline Scientific Analysis (OSA) 
software distributed by ISDC (' Courvoisier et al.l l2003l) is op- 
timized for point sources, and is not well suited for slightly 
extended sources. We present here a method based on Pixel 
Illumination Fraction (PIF) which extracts the properties of 
slightly extended sou rces with a coded mask instrument (see also 
iRenaud et al.]|2006bl) . 

For a coded mask instrument, the sky images and, in partic- 
ular the mosaic image studied in the previous section, are pro- 
duced by the backprojection of the shadow patterns cast by the 
sources in the field of view on the plane of the sky. The shadow 
pattern produced when observing a FOV containing n sources 
is a superposition of the shadow patterns of all the individual 
sources. 



(6) 



1=1 



where fj and PIF,(jic, y) are the flux, respectively the shadow pat- 
tern (called "Pixel Illumination Fraction") of the /th source and 
B{x, y) is the background in the pixel with coordinates (x, y). The 
Pixel Illumination Fraction gives the fraction of each pixel of the 
detector illuminated by the source. For a pixel that is completely 
in the shadow of the mask, the PIF will be 0, whereas in the case 
of a fully illuminated pixel, the PIF will be equal to 1 . It is un- 
derstandable that the PIF of an extended source is different from 
that of a point source, since some pixels might be illuminated 
by only a fraction of the total extension of the source, which 
cannot happen for a point source. Thus to describe an extended 
source properly, one has to create an appropriate model for it. 
Our method to create such a model is the following: we create a 
grid of positions on the sky covering the extended source, com- 
pute the PIF for all the positions of this grid, and then average 
the PIFs, weighted by a model surface brightness for the source, 
e.g. a spherical isothermal yS-profile, 



/(r) 



1 



Fig. |5] shows the PIF for an on-axis source, in the case of a 
point source and for an extended source described by the model 
ofEq.|7]wiflia = 30'. 




(7) 



Fig. 5. The Pixel Illumination Fraction (PIF) of the IBIS/ISGRI 
instrument for a source in the center of the field of view 
(white=l, black=0): left: for a point source; right: for an ex- 
tended source with a surface brightness described by an isother- 
mal y6-profile, with a source size a = 30 arcmin. 



Since the standard OSA tools always use the PIF for a point 
source, in the case of an extended source they do not extract 
fluxes and spectra correctly. To extract the flux of an extended 
source, we created a tool that fits the detector image (or shadow- 
gram) with any kind of PIF, either point-like or extended with 
a given surface brightness model. Our tool creates a correct PIF 
for every source in the FOV, and then fits the shadowgram to 
the model described in Eq. (|6]l, where 1 or more sources might 
be extended. To check the coherence of our method, we have 
simulated the shadowgram one can expect with 2 sources, one 
extended in the middle of the FOV with a flux Fx - AQ counts 
per fully illuminated pixel (cpp hereafter) and one point-like in 
the corner of the FOV with a flux F2 - 80 cpp, with a gaussian 
background, and extracted the fluxes of the 2 sources with our 
tool. The results gave us Fi = 40.0 + 0.7 cpp and F2 = 80.4 + 1 .2 
cpp, which shows that our method is indeed extracting fluxes 
properly. 

This method can also be used to extract spectra: we analyzed 
the data with the standard OSA tools to create shadowgrams in 
all the energy bands desired for the spectrum, and then used our 
fitting tool to extract a flux with the correct PIF in all the energy 
bands, and reconstruct a spectrum (see Sect.|5]l. 

It is important to note that the use of a PIF-based method 
for weak extended sources could be complicated because of the 
specifics of operation of the INTEGRAL satellite. Namely, each 
INTEGRAL observation is split into several kilosecond-long 
intervals of continuous data taking, called Science Windows 
(ScWs). The statistics of the signal from a weak source in each 
ScW is below the background statistics. Moreover, the statisti- 
cal properties of the low-statistic signal in the ISGRI imager are 
not well known. This means that the use of gaussian statistics 
for the fitting of the shadowgrams is not completely correct, and 
can give wrong results. On the contrary, the method of analy- 
sis based on the mosaic image is better suited for the analysis 
of weak extended sources. In any case, the two methods should 
be considered complementary, and for this reason we present re- 
sults based on both methods in this paper. 
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Fig. 6. Substraction between the mosaic image in the 18-23 keV 
band and the model obtained with the fit of the data with the sum 
of NGC 4874, EXO 1256+281 and NGC 4889. While the center 
of the emission is well described, the outer parts of the image 
show a deviation of up to 5.5cr from the model. 



4.2. Hypothesis of multiple point sources 

One possible explanation for the shape of the mosaic image (Fig. 
[U is that one or several point sources contribute to the observed 
flux. We analyse this possibility here. Candidate sources are the 
cluster itself and bright AGNs at close angular distances. The 
three brightest AGNs in soft X-rays are the central radio galaxy 
NGC 4874, NGC 4889 and the QSO EXO 1256+281. Fig. IT] 
shows the position of these 3 point sources on the ISGRl image. 

We used the method described in Sect. 14.11 to fit the shad- 
owgram in the 18-23 keV band by the PIF of 3 point sources at 
the position of the sources given above. We extracted the flux of 
each of these 3 sources, and created a model of the deconvolved 
image with the flux found. 

Fig.|6]shows the result of the substraction of the fitted model 
from the mosaic image in the 18-23 keV. We can see on the im- 
age that the outer parts of the source are not well described by 
the model, specially the region at the bottom right and the left, 
where 5.5cr excesses are observed in the mosaic compared to 
the model image. This analysis suggests that the emission seen 
by INTEGRAL is more extended than what can be explained by 
the superposition of the brightest AGNs in the cluster. 

To confirm this, we fit the data with 2 fake point sources 
placed along the major axis of the ellipse, and computed the de- 
tection significance for different distances of the 2 fake sources. 
In the case of a single point source, we would expect the detec- 
tion significance to drop when we increase the distance between 
the 2 fake sources, whereas in the case of 2 sources whose PSF 
overlap, we expect the detection significance to peak at the dis- 
tance between those 2 sources. The position of the 2 sources for 
which the detection significance peaks will therefore allow us 
to compare this result with the position of known point sources 
inside the cluster. If there is a possible point source counterpart 
in soft X-rays, this will give an indication that contamination of 
the hard X-ray flux by point sources is likely. The result of this 
computation is shown in Fig.|7] 

We can see on Fig. [T] that the detection significance peaks 
at a distance of 17 ± 2 [arcmin], which is fully compatible with 
the result obtained from the mosaic in Sect. 13. II Since we expect 
the thermal emission from the cluster to dominate in this energy 



13 - 



12 - 



5 10 15 20 25 30 35 

Distance between points [arcmi n] 

Fig. 7. Computation of the total detection significance in the 18- 
23 keV band as a function of the distance between 2 fake sources 
along the major axis of the ellipse. 




5.5 - 



5 10 15 20 25 30 35 

Major axis of the ellipse [arcmin] 

Fig. 8. Detection significance of the source computed with an 
ellipse-shaped PIF in the 18-23 keV band as a function of the 
major axis of the ellipse. 

range, we would expect the two points to be located on the centre 
of the cluster (that lies 1' south of NGC 4874) and one of the 
other AGNs. However, the distance between the centre of the 
cluster and EXO 1256+281 (4.6') or NGC 4889 (8.9') is at least 
a factor of 2 smaller than the distance we found. This is a strong 
argument in favour of extended emission. 

4.3. Extended emission analysis 

We used the same method to analyse the ISGRI data with the as- 
sumption that the emission is indeed diffuse emission. We used 
the angle of the ellipse found from the mosaic analysis (see Sect. 
13.11 ) to create PIFs of an extended source with the method ex- 
plained in Sect. 14.11 weighted by a surface brightness given by 
Model 3 (see Eq. (|3]l). We used our tool to fit all the data with 
this model for several different values of cri and 0-2 with a ratio 
^ - 1.43 fixed by the results of the image fitting, and finally 
computed the detection significance of the source for all the dif- 
ferent models. The result is shown in Fig. [8] 

We can see in this Fig. that the detection significance also 
peaks at a major axis size of about 17 [arcmin], which is again 
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comparable to the results obtained from the mosaic. We can see 
that the maximum detection significance is lower than in the case 
of 2 point sources, but this does not mean that the model of 2 
point sources gives a better representation of the data: indeed, 
the eiTor n the fit increases with the size of the source, because 
the contours of the smallest holes of the mask become more and 
more unclear (see Fig.|5]l, and thus the imaging method is less ac- 
curate. It is thus irrelevant to compare directly Figs. [T] and [8] but 
both show independantly the fact that the source is not point-like, 
and provide a measurement of the apparent size of the source in 
hard X-rays. 

From this analysis, we conclude that if we use a source size 
smaller than the limit of 17', we lose a part of the total flux of 
the cluster, and thus the detection significance increases with the 
size. In the opposite case, we collect more background, so the 
detection significance starts to drop. The best estimate of the to- 
tal flux of the cluster is therefore given at a major axis size of 
17 [arcmin], which gives F,„, - 0.31 + 0.04 counts/s in the 18- 
30 keV band. From now on, we will use this model to extract 
fluxes and spectra. 

We also performed another complementary analysis: we fit- 
ted the data with 4 sources, i.e. the extended model described 
above and the 3 AGNs. The fit converges to a solution that puts 
80% of the flux in the extended emission. This model cannot 
be used to extract a flux, because the fluxes of the 4 sources 
become strongly anti-correlated. However, it means that the ex- 
tended source model describes the data better than the sum of 
point sources. 

5. Spectral analysis 

To extract the spectrum of the source correctly, we used the 
method described in Sect. 14.11 We used the hard X-ray shape 
of the source extracted from the INTEGRAL image to create a 
PIF covering the whole size of the source, and extracted the 
flux from the shadowgram of each pointing in 3 different en- 
ergy bands: 18-23, 23-30 and 30-40 keV (the source is not de- 
tected at higher energies). We then performed a weighted sum 
over the flux extracted from all pointings to get a total spec- 
trum in these 3 energy bands. We also extracted the XMM/PN 
spectrum of the cluster with the background substraction method 
from the Birmingham group, in a region chosen such that the 
XMM-Newton/INTEGRAL intercalibration factor is equal to 1 . 
Finally we fitted this spectrii m in XSPEC with the MEKAL 
model dKaastra & Mewel2000l) . We extrapolated the fitted model 
to higher energies and compared it to the data points obtained 
with ISGRI. The result is shown in Fig.|9] 

Because of the very low statistics at energies above 30 keV, 
we are not able to confirm or deny the presence of a non-thermal 
hard X-ray excess emission. Indeed, the extended nature of the 
source makes it difficult to extract a significant spectrum up to 
high energies, since the already low statistics is spread over sev- 
eral sky pixels. A longer exposure time is therefore required to 
make conclusions on the presence or not of a hard X-ray excess 
emission above 30 keV. 

Although there is strong evidence that EXO 1256-1-281 can- 
not explain the INTEGRAL SW extension, the angular distance 
between the fitted position and EXO 1256-H281 (6.1') is close 
to the half-width of the ISGRI PSF Hence, we further investi- 
gate this point by extracting the spectral properties of this ob- 
ject, assuming the identification of the additional point source 
with EXO 1256+281. 

To this end, we fitted the excess in the ISGRI image by a 
point source at the position of EXO 1256+281, and extracted 



Combined XMM/PN and INTEGRAL/ISGRI spectrum 
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Fig. 9. Combined XMM/PN and INTEGRAL/lSGRl spectrum. 
The PN spectrum is fitted by a MEKAL model at a tempera- 
ture kT - 7.9 + 0.1 keV. The bottom plot shows the residuals of 
the best fit compared to the data. There is no obvious deviation 
to this spectrum from the ISGRI data points up to 40 keV. 
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Fig. 10. XMM-Newton spectrum of the quasar EXO 1256+281 
with one ISGRI data point obtained by fitting the south-west ex- 
cess in hard X-rays by a point source at the position of EXO 
1256+281. The ISGRI upper Hmit in the 30-50 keV band is 
shown. The spectrum is fitted by the sum of a power law (dashed 
red) and a heavily absorbed component dominating the flux at 
higher energies (dashed green). 



the flux for this source. Fig. [TOl shows the XMM-Newton spec- 
trum of EXO 1256+281 with the ISGRI point extracted using 
this method. The XMM-Newton spectrum is well fitted by a sim- 
ple powerlaw with the photon index F = 1.63 + 0.07 and 0.5- 
10 keV flux 5 X 10"'^ ergs cm"- s"'. The flux of the source in 
the INTEGRAL energy band is much higher than the extrapola- 
tion of the powerlaw found in XMM-Newton data. This implies 
the presence of an additional spectral component which domi- 
nates the source above 10 keV. To model this component, we 
added a heavily absorbed power law typical of Seyfert 2 galax- 
ies to the fit. Taking into account the upper limit in the 30-50 
keV band, the fit gives = (4.0 + 1.7) x lO^"* cm^^, and F > 3.0 
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to match the upper limit. This index is much steeper than the 
unified Sey fert 2 spectral index (F - 1.79 with a dispersion 
cr = 0.23, iRisali ti (2002)), and hence, the properties of this 
source would be very unusual for a Seyfert 2 galaxy. From this 
statement together with the imaging arguments presented in the 
previous sections, we conclude that the contribution of known 
point sources to the observed flux is very unlikely between 10 
and 40 keV. 

Assuming that the SW excess is due to the presence of dif- 
fuse emission, which gives the best representation of the data, 
one can try to constrain the properties of the gas needed to ex- 
plain the emission from this region. We attempted to make a joint 
fit of the 1-10 keV spectrum extracted from a circle of a radius of 
6' centered at the position of the SW excess and the 18-50 keV 
spectrum extracted from the SW region from the INTEGRAL 
mosaic image in the same way as in the case of an additional 
point source. Assuming that all the flux from the SW region 
comes from a higher temperature plasma, we fitted the data with 
the thermal bremsstrahlung model. This results in a temperature 
of kT = 12 + 2 keV. The corresponding estimate of the emission 
measure (EM) is 0.16 < EM < 0.26 cm"^ pc, which is rea- 
sonable for external regions of the cluster. We can thus conclude 
that the presence of a hotter region (10 keV< kT < 14 keV) can 
explain the extension found in the ISGRI mosaic image. 



6. Discussion 

In this paper, we have used the IBIS/ISGRI instrument on board 
the INTEGRAL satellite to investigate the hard X-ray emission 
from the Coma cluster. We presented a method based on Pixel 
Illumination Fraction (PIE) to analyse extended sources with a 
coded mask instrument (section |4n ), and we have shown that the 
Coma cluster indeed appears like an extended source for ISGRI. 
Assuming that the emission seen by INTEGRAL is extended, we 
have compared the ISGRI mosaic image with the soft X-ray im- 
age from XMM-Newton, and shown that there is a displacement 
between them: the INTEGRAL image is displaced towards the 
south-west, i.e. in the direction of the NGC 4839 group, that is 
cuiTently merging with the main cluster. 

The origin of the extended emission from this region is not 
clear, but we have investigated two possible explanations for the 
excess in the image: an additional heavily absorbed point source 
embedded in the cluster, and an extended region where large- 
scale shocks occur. 

The first model for the hard X-ray excess in the South- West 
region of the cluster is the presence of a highly absorbed addi- 
tional point source appearing at higher energies. We have used 
the imaging capabilities of the instrument to investigate this pos- 
sibility: we have shown that the shape of the residual image after 
substraction of the XMM-Newton surface brightness profile from 
the ISGRI image does not coincide with any known X-ray point 
source, and that the only possible candidate EXO 1256+281 is 
located more than 6 arcmin away from the best position found in 
the INTEGRAL image, which makes it an unlikely counterpart. 
We have also extracted the soft X-ray spectrum of this source 
and the flux of the south-west region in INTEGRAL data, and 
shown that this spectrum is not compatible with a highly ab- 
sorbed Seyfert II galaxy. Indeed, the source is not detected in 
the 30-50 keV band, which is the most sensitive energy band 
of ISGRI. This implies a spectral index F > 3.0, which is too 
steep for a highly absorbed Seyfert 2 galaxy. As a conclusion, 
we claim that contribution of a very hard point source embedded 
in the cluster to the observed spectrum is highly unlikely. 



We cannot exclude the possibility that the South- West excess 
in hard X-rays is due to one or few unknown sources that would 
emit predominantly in this energy band. If this is the case, these 
sources must be highly absorbed {hh > 4 x 10^^ cm^), and have 
a steep spectral index F > 3.0. We note that these characteris- 
tics are unlikely, because most of the highly absorbed sources 
discovered by INTEGRAL show a much harder spectrum. 

In the scenario of a merging event between the Coma cluster 
and the NGC 4839 group, we expect a shock front to be cre- 
ated in the region where the gas of the two clusters collides. Our 
imaging analysis shows that the hard X-ray emission seen by 
INTEGRAL is extended in the direction of the NGC 4839 group, 
which is a good indication that the emission we see is indeed 
coming from a region where large-scale shocks occur If this ex- 
planation is correct, we expect the plasma in this region to be 
hotter, and hence to have a h arder spectrum in X-ra ys. The tem- 
perature map of the cluster (Ne umann et al.l 120 03) shows a re- 
gion that has low surface brightness in X-rays, but is the hottest 
region of the cluster (kT > 10 keV). This region coincides with 
the position of the hard X-ray South- West extension found by 
INTEGRAL, and we can thus associate the hard X-ray excess in 
this region discussed in Sect.|5]with emission from a very hot re- 
gion of the cluster {kT < 14 keV). This result is consistent with 
the temperature fo und in the merging re gion of the distant cluster 
CI J0152.7-1357 (Mau ghan et al.l20()3i') . which shows that such a 
high temperature is possible and might indeed be the signature of 
a merger Assuming that the hot region is roughly spherically sy- 
metric and has an angular size of ~ 6' at a distance of 100 Mpc, 
we deduce that the density of the hot gas is n/,,,, ~ 10"^ cm"^. 
Given that the South- West excess in the ISGRI image appears to 
be extended, INTEGRAL data appear to confirm this scenario. 
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